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Evidence for local charge fluctuations linked to a charge disproportionation of the Bi ions in the 
distorted lattice of superconducting LaOi- x F 2 ,BiS 2 is presented. In-plane short-range distortions 
of sulfur atoms up to 0.3 A in magnitude break site symmetry and create two distinct environ¬ 
ments around Bi. Out-of-plane motion of apical sulfur brings it closer to the La-O/F doping layer 
with increasing x that may lead to a charge transfer conduit between the doping layers and the 
superconducting BiS 2 planes. The mechanism for superconductivity may arise from the interplay 
between charge density fluctuations and an enhanced spin-orbit coupling suggested theoretically, 
that induces spin polarization. 

PACS numbers: 61.05.F-, 74.62.Dh, 74.70.Ad, 75.50.Bb 


The rich physics of the Bismuth Sulfide (BiS 2 ) based 
superconductors includes a possible charge density 
wave (CDW) instability@, 9$, atomic modes leading to 
phonon softeningQ and spin-orbit (SO) interactions with 
a hidden spin polarization [10]. Much of the recent in¬ 
terest in superconductors has been on materials whose 
parent ground state is magnetic i.e. the Fe and Cu based 
systems. The BiS 2 superconductors are not magnetic, 
possibly of the conventional type but with strong elec¬ 
tronic correlations Em. Density functional theory cal¬ 
culations suggested that spin polarization is present even 
though the crystal lattice is centrosymmetric @[k| , aris¬ 
ing from coupling of the electron orbital to its spin. Such 
SO coupling has been at the heart of spintronic physics, 
exemplified at the interfaces of films like LaAlOa/SrTiOa 
EH: where a large Rashba SO couplin g is important in 
stabilizing the superconducting state [l3|. Given its 
quasi-two dimensional crystal nature, BiS 2 may exhibit 
SO coupling through the Dresselhaus effect [10|, the 3- 
dimensional analogue of the Rashba effect. This can 
happen without breaking crystal inversion symmetry but 
rather be driven by atomic site asymmetry El- 

In the two families of BiS 2 superconductors, 
Bi-tO^SO^i-a, fi] and LnOi_ s F a; BiS 2 (Ln = La, Nd, Pr, 
Ce, and Yb) [l4, [l5[ discovered thus far, vacancies and 
crystal defects are common. While the SO 4 layer carries 
vacancies in the former compound, the LnOi_ 3 ,F a; BiS 2 
system shows a peculiar relation of its crystal structure 
to superconductivity. The parent compound shows high 
crystallinity but upon doping, T q is maximized when 
crystal disorder is induced by high pressure annealing. 
The relation of disorder and superconductivity has been 
intensely studied in cuprates [16| and, earlier, in amor¬ 
phous solids Emm- Disorder shortens the coherence 
length, enhances thermal fluctuations and induces vortex 
pinning but without long-range ordering into a vortex lat¬ 
tice [lg, |19( . Such a state has been argued to correspond 


to a vortex glass [20|. However, the structural disorder 
observed in LnOi-^F^B^ is different. 

In this system, the dominant carriers are from the Bi 
6p orbitals that strongly hybridize with the S 3p near the 
Fermi surface. A CDW instability arises from a ( 7 r, 7 r, 0 ) 
mode which corresponds to in-plane displacements of S 
atoms (0.18 A) around the M point[ 9 ]. At the same 
time, SO coupling has a considerable effect on the elec¬ 
tronic band structure around the X point of the Brillouin 
zone. Without SO coupling, four bands cross below 
in the optimally doped composition, x ~ 0.5. With SO 
coupling, only two bands cross below E^. SO coupling 
originates from the inversion asymmetry and creates an 
effective magnetic field that lifts the spin degeneracy [Hj. 
Are the new Bi-based superconductors conventional BCS 
type or a new class of non-BCS type, yielding a new route 
to unconventional superconductivity? Here, we provide 
evidence for the latter by analyzing neutron scattering 
data. Our results provide the experimental evidence 
for in-plane short-range S distortions (~ 0.3 A in magni¬ 
tude) that lead to charge fluctuations around Bi. Crystal 
strain develops along the compressed c-axis but the out- 
of-plane S atoms move closer to the doping layers. The 
local distortions present in this system may be an exper¬ 
imental realization of a vortex pinning mechanism in the 
glassy state of a superconductor. 

The quasi two-dimensional layered crystal structure of 
LaOBiS 2 (see ref [ 3 - 5 ] for sample preparation) shown 
in Fig. 1(a) consists of the superconducting BiS 2 bilay¬ 
ers sandwiched between insulating La 2 C >2 blocking layers. 
The Bi ions form a square pyramidal edge shared lattice, 
coordinated by four in-plane (SI) and one out-of-plane 
(S2) sulfur ions. The plane is buckled with a Sl-Bi-Sl 
bond angle of ~172.1°. In the tetragonal symmetry, the 
base of the pyramid is a perfect square, and the Bi-Sl 
bonds are of the same length, 2.876 A. The out-of-plane 
Bi-S 2 bond length of 2.476 A is significantly shorter. The 
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substitution of oxygen (O 2- ) with fluorine (F _ ) in the 
blocking layers introduces electron carriers, and while the 
crystal symmetry, PJ^/nmm, does not change, the c-axis 
contracts by 4 % and the a-axis expands by less than 0.2 
%. The samples become superconducting under high 
pressure annealing. The BiS 2 layer buckling is reversed 
with doping, in which the Bi ion moves away from the 
S2 atom, and the Bi-S2 bond length increases to 2.65 A 
by x = 0.5. The Sl-Bi-Sl bond angle becomes 184.8° 
with a Bi-Sl bond of 2.87 A. This contradicts theoretical 
predictions that expect the lattice to become flat with 
doping [2^ . 

The centrosvmmetric space group does not change 
with doping |2|, even though significant Bragg peak 
broadening along with diffuse scattering are observed in 
the neutron diffraction pattern (see Fig. 1(b) which is a 
comparison of the x = 0.0 and x = 0.5). It is perplexing 
how such a distorted crystal lattice is conducive to su¬ 
perconductivity. Earlier studies showed that the broad¬ 
ening affects Bragg peaks with nonzero /-indices while 
Bragg peaks with l = 0 are much sharper as reported 
in Ref. j23|. Fig. 1(d) is a comparison of two nuclear 
Bragg peaks, the (110) and (114). The (114) Bragg peak 
and others with nonzero l become quite broad in the 
high pressure treated samples. This has been attributed 
to c-axis strain [24], arising from stacking faults in the 


presence of partial dislocations from imperfect stacking 
of crystal planes. 

To investigate the origin of the broadening and dif¬ 
fuse scattering, here we use a real-space approach to ex¬ 
tract the atom specific characteristics of the distortions 
and their real-space configuration. Bragg peak broad¬ 
ening and diffuse scattering can arise from other form 
of defects besides stacking faults: the size mismatch be¬ 
tween O and F can contribute to Huang scattering, for 
instance, while thermal diffuse scattering (TDS) due to 
phonons and static displacements of ions can addition¬ 
ally contribute to the diffuse scattering. We previously 
estimated the Debye temperature to be around 200 K for 
the x = 0.5 [25l | which indicates that TDS may not be 
insignificant in the temperatures of interest here, < 10 
K. The real-space atomic structure obtained by Fourier 
transforming the total structure factor that includes both 
Bragg and diffuse scattering is represented in the pair 
density function (PDF), p(r). The p (r) is plotted in 
Fig. 2(a) for four compositions. This analysis provides 
the local arrangement of atoms and is very sensitive to 
short-range distortions. Further details can be found in 
Refs: 


I27| . The p ( r ) consists of correlation peaks 
that correspond to the probability of finding a particular 
pair of atoms, as labelled in the figure, at a given distance 
in space. 

While in reciprocal space Bragg peak broadening is 
very pronounced as seen in Fig. 1(c), in real space, the 
PDF correlation peaks change shape with x , an unusual 
effect that indicates a local rearrangement is taking place 


that may be driven sterically and/or energetically. The 
crystal symmetry is invariant under doping and the dif¬ 
ferences observed in the PDFs of LaOi_ a; F a ,BiS 2 cannot 
be explained by the substitution of F _ for 0 2 ~ as the two 
have very similar neutron scattering lengths and nomi¬ 
nal ionic sizes. Starting with the parent compound, the 
PDF data is compared to a model p(r) mo d, calculated 
using the unit cell dimensions and atomic coordinates of 
the average structure. In Fig. 2(b), p(r) exP (symbols) 
is compared to the model p(r ) moc j (dash blue line - av¬ 
erage). The overall agreement between the two curves is 
good, however very clear differences are observed espe¬ 
cially around 2.5 -3.5 A that are not due to systematic 
or statistical errors. Some ripples can be seen around 
the first PDF peak due to termination effects resulting 
from the truncation of the Fourier transform at a finite 
momentum transfer but the aforementioned differences 
are well above the error. 

In-plane SI displacements can reproduce the local con¬ 
figuration that best fits the p(r) exp of Fig. 2(b). The 
displacement modes are shown in Figs. 2(c) and 2(d). 
In 2(c), the displacement mode is of the breathing type. 
Only SI is displaced while all other atom coordinates 
are kept the same as in the average structure. A new 
p(?')mod is calculated using the breathing mode model 
and is compared to the p(r) exp also in Fig. 2(b) (solid 
red line - local). In spite of its simplicity, this model 
fits the data quite well, especially in the 2.5 to 3.5 A 
range. The breathing mode shown in Fig. 2(c) creates 
short and long Bi-Sl bonds that affect the charge distri¬ 
bution around Bi. The distortion creates a double-well 
displacement potential of the Bi-Sl bonds. The data can 
be reproduced equally well using the displacement mode 
of Fig. 2(d). In this model, the SI atoms are displaced 
in a ferrodistortive type mode, and create short and long 
Bi-Sl bonds in an asymmetric way. Both modes are 
equally likely as they give rise to identical magnitude of 
bonds. It is conceivable that both the ferrodistortive and 
breathing type modes are locally present. The asymmet¬ 
ric coordination is reminiscent of the environment of Bi 
in Cu 4 Bi 5 Sio @ due to the presence of a lone pair. 

In the superconducting state, similar local modes are 
likely. The p(r) exp of x = 0.5 is shown in Fig. 3(a) 
and is compared to the p(r) moc j for the average model of 
x = 0.5 (dash blue line). Differences are observed be¬ 
tween the two in the same region of space as in the parent 
compound even though the shape of the PDF peaks are 
different. Just as in the parent compound, the breath¬ 
ing and ferrodistortive modes can fit the data quite well 
(p( r )mod’s shown in solid lines) but the magnitude of dis¬ 
tortion is different. Shown in the inset of Fig. 3(a) are 
the Bi-Sl and Bi-S2 partial PDFs obtained from the local 
and average models for x = 0 and x = 0.5. It can clearly 
be seen that the average Bi-Sl bond of 2.8 A is split to 
two locally in both the x = 0 and 0.5 while the split is 
larger in the x = 0 than in x = 0.5 sample. A third small 
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peak is observed in close proximity to this bond that is 
from Bi-Sl bond lengths across planes. 

The charge fluctuation inferred from the long and short 
Bi-Sl bond lengths is present in all compositions leading 
to the x = 0.50 which is the only superconducting com¬ 
position in our series. With doping, the Bi - SI bonds 
differ in length as shown in the schematic of the Bi tetra- 
hedra of Fig. 3(b). We can also see that upon entering 
the superconducting phase, the apical S2 atom gets fur¬ 
ther away from the planes and closer to the La(0/F) 
charge layers. The S2 atom can mediate the transfer 
of charge from the doping layers to the conduction BiS 2 
layers just like the chain oxygen does in YBa 2 Cu 307 _a. 
With temperature, the split of the Bi-Sl bonds in the x 
= 0.5 changes little except on approaching T<y where the 
split becomes larger (see Fig. 4(a)). On the other hand, 
it is clear that with cooling, the height of the S2 atom 
gets closer to the donor planes (Fig. 4(b)). This is also 
consistent with the temperature dependence of S2 in the 
average structure. 

The long and short Bi-Sl in-plane bonds indicate a 
variance in the Bi valence due to charge disproportion¬ 
ation. The proposed charge fluctuations are shown in 
the schematic of Fig. 3(c) for the BiS 2 planes. Two dis¬ 
placement modes of the breathing type are tested: given 
that there are two planes at 2 ~ 1/3 and 2 ~ 2/3, in 
one mode the SI distortions are in the (x, -x) direction 
in each plane respectively, while in the other, the distor¬ 
tions are in the (x, -y) direction in each plane. Both 
modes yield the same local bond order but might be dif¬ 
ferent energetically. Note that within one unit cell, the 
two BiS 2 planes are rotated by 90 degrees. The (x, -x) 
only breaks the 4-fold rotational symmetry to a 2-fold 
while the inversion symmetry remains. The (x, -y) dis¬ 
tortions break the 4-fold and 2-fold rotational as well as 
the inversion symmetries. 

The superconducting properties depend on structural 
defects and on their arrangement, and knowing the av¬ 
erage or ideal structure is not sufficient [30j. What can 
the local structure tell us about the crystal state of this 
new superconductor? It is clear that superconductiv¬ 
ity is enhanced with the quenched disorder in this sys¬ 
tem and SI and S2 distortions are important. Theo¬ 
retical works have shown that a vortex-glass supercon¬ 
ductor may be present in the bulk of disorder supercon¬ 
ductors. In the presence of disorder, superconducting 
vortices are pinned, preventing the formation of vortex 
lattices 0- Thus even though the BiS 2 planes are dis¬ 
torted in the non-superconducting and superconducting 
samples, long-range order exists in the former as evi¬ 
denced from the Bragg structure, and in the absence 
of pinning centers, the system can be in a vortex fluid 
phase with strong fluctuations preventing pairing. On 
the other hand, in the superconducting sample, the dis¬ 
order is strong not only due to the distortions of the 
BiS 2 planes but also due to the c-axis strain and stack¬ 


ing faults, allowing the system to enter the vortex glass 
thermodynamic state. The stripes of distortions pro¬ 
posed in Fig. 3(c) are short-range and the presence of 
stacking faults break the periodicity thus preventing any 
long-range order. It is possible that such modes are the 
cause of vortex pinning. 

The doping of F _ for O 2 ^ changes the oxidation state 
of Bi nominally from 3+ to 2+ while introducing elec¬ 
tron carriers. A simple estimate of the Bi 3+ - S 2- bond 
length in the x = 0 yields a value of ~ 2.8 A that is close 
to the bond length determined from the average struc¬ 
ture. However, locally, the Bi 3+ - S~ 2 bond lengths are 
split, as shown above. We cannot distinguish whether 
this is the result of a breathing or a ferrodistortive type 
mode because they yield the same bond magnitude. The 
leading instability mode consists of displacements of SI 
in the superconducting plane and out-of-plane S2 distor¬ 
tions. The in-plane distortions get smaller in the high 
pressure annealed x = 0.5 superconducting samples while 
the out-of plane z-motion of S2 increases. The two pro¬ 
posed modes will both create unequal charge distribu¬ 
tions. However, such charge fluctuations are local and 
cannot produce a charge density wave. Even though this 
system is not magnetic, theoretical calculations showed 
that SO coupling leads to a hidden spin polarization that 
has a considerable effect on the electronic band structure 
around the X-point of the Brillouin zone. How the insta¬ 
bility mode proposed here affects spin polarization will 
be interesting to investigate. In general, superconductiv¬ 
ity in the presence of strong structural disorder might be 
a playing field for discovering exotic states. 
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Fig. 1. (a) The crystal structure of LaOi_ x Fj.BiS 2 
with the BiS tetrahedra. (b) The neutron powder diffrac¬ 
tion pattern (black symbols) of LaOBiS 2 at 6 K was col¬ 
lected using NOMAD at the Spallation Neutron Source 
of Oak Ridge National Laboratory. The red, green and 
blue solid lines represent the calculated intensity, back¬ 
ground and the difference between the observed and cal¬ 
culated intensities, respectively, (c) The neutron powder 
diffraction pattern of LaOo. 5 Fo. 5 BiS 2 at 2 K has its in¬ 
tensity considerably reduced in comparison to the parent 
compound of (b). At the same time, the background 
(green solid line) has increased. (d) The comparison 
of the (110) and (114) Bragg peaks between the x = 0 
and x = 0.5 compositions shows significant changes of 
the non-zero l peaks with x. 


Fig. 2. (a) A plot of the composition dependence of the 
data corresponding to the local structures of x = 0, 0.2, 
0.3 and 0.5. (b) The PDF of the local structure for the 
parent compound at 6 K (black symbols) is compared to 
the average (dashed blue line) and local (solid red line) 
model PDFs. The agreement factor (A/ octor ) between 
the local model and data is 0.156 while between the av¬ 
erage model and data it is 0.283. (c) A schematic of 

the BiS 2 plane showing the breathing type mode, (d) A 
schematic of the BiS 2 plane showing the ferro-distortive 
type mode. The displacements of SI are either in the x- 
or y-direction. 

Fig. 3. (a) The PDF of the local structure for x 

= 0.5 at 2 K (black symbols) is compared to the aver¬ 
age (dashed blue line) and local (solid red line) model 
PDFs. The Af actor between the local model and data 
is 0.18 while between the average model and data it is 
0.239. The inset is a plot of the partial PDFs for Bi-Sl 
and Bi-S2 for x = 0 (top) and x = 0.5 (bottom). The 
partial of Bi-Sl for both the average and local models is 
shown while the partial of Bi-S2 is only shown for the lo¬ 
cal model, (b) A schematic of the Bi-S tetrahedron with 
the bond lengths determined from the refinement of the 
local structures of 0 < x < 0.5. (c) Stripes of charge fluc¬ 
tuations in the two BiS 2 planes of the crystal structure. 
The (x,-x) and (x,-y) refer to the SI coordinates in the z 
~ 1/3 and z ~ 2/3 planes, respectively. 

Fig. 4. (a) The temperature dependence of the 

short (square) and long (circle) local Bi-Sl bond lengths 
is contrasted to the Bi-Sl average bond length (triangle), 
(b) The temperature dependence of the S2 atom from the 
oxygen plane obtained from fitting the data and from the 
Rietveld refinement of the average structure. 
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(c)LaO 05 F 05 BiS 2 , T=2K 
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Fig. 1, Athauda et al. 
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Fig. 2, Athauda et al 
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Fig. 3, Athauda et al. 
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